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Abstract 

A neutral network is a subgraph of a Hamming graph, and its principal eigenvalue determines its robustness: the 
ability of a population evolving on it to withstand errors. Here we consider the most robust small neutral networks: the 
graphs that interpolate pointwise between hypercube graphs of consecutive dimension (the point, line, line and point 
in the square, square, square and point in the cube, and so on). We prove that the principal eigenvalue of the adjacency 
matrix of these graphs is bounded by the logarithm of the number of vertices, and we conjecture an analogous result 
for Hamming graphs of alphabet size greater than two. 
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The eigenvalues of neutral networks—subgraphs of Hamming graphs—is a fascinating subject, yet one which seems to 
have received little attention from the mathematics community. A recent surge of scientific interest has been motivated 
by advances in the theory of neutral evolution DM in which the evolution of a mutating population is captured by 
spectral properties of its underlying neutral network 0. 

A genome is the set of all genotypes, or a-ary strings, of length d and alphabet size a. Typically a is small: a = 2 
(hydrophilic and hydrophobic), a — 4 (nucleic acids) or a - 20 (amino acids). On the other hand, d can range 
from 3 (codons) to 10 s (chromosomes). We represent the genome by a rZ-dimensional Hamming graph ll t j „ = (K a ) d = 
(K„U ... OK a ), where K a is the complete graph on a vertices and □ is the Cartesian product (4). Each of the a d vertices 
in the Hamming graph corresponds to a genotype, and two vertices share an edge if the genotypes differ by a single 
mutation (Hamming distance one). A neutral network is the set of genotypes with the same phenotype (observable 
characteristics); it is just a subgraph of //,/„. In this Note we use neutral network and phenotype interchangably. 

Assuming the phenotype has achieved high fitness, adjoining phenotypes will have a relatively negligible growth 
rate and act as effective absorbing boundaries. Now consider a mutating population on the neutral network. That 
portion which mutates off of it will be lost, whereas that portion which stays on will survive. The robustness r of 
a neutral network is the long-term probability that a randomly selected individual mutating in a random direction 
survives. It is the principal eigenvalue A of the adjacency matrix of the neutral network divided by the number of 
directions for mutation: r = A/(d(a - 1)). This Note is about the maximum robustness that a neutral network can have, 
which for small neutral networks are themselves Hamming graphs and interpolations between them, shown for a = 2 
in Figure 1. Surprisingly, this ceases to be true for larger neutral networks on longer sequences (d = 19 and above), 
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Figure 1: The first 16 “bricklayer’s graphs,” which are the Hamming graphs and the interpolations between them. Below each is the principal 
eigenvalue of its adjacency matrix. 
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which we demonstrate by way of counterexamples in the form of star graphs. 

This Note is divided into six parts. In the first part, we derive the relation between the robustness r of a neutral 
network and its principal eigenvalue A\ r - A/(d(a - 1)). In the second part, we present our main result, A„ < log 2 n, 
where A„ is the eigenvalue of the nth bricklayer’s graph (see Figure 1). We prove it for all values of n apart from 
n — 2 d ± 1 using a geometric staircase argument. In the third part, we set the stage for the remainder of the proof by 
constructing polynomials with eigenvalues A 2 d±\ as roots. In the fourth and fifth parts, we bound the eigenvalues A 2 «_ \ 
and A 2 d+i . In the sixth part we conjecture a generalisation of our main result for higher a and conclude with some 
extensions. 

A secondary motivation of this paper is to introduce mathematicians to the connection between neutral networks 
and spectral graph theory and to encourage them to extend this work. 


1. Relation between the robustness and the principal eigenvalue 


In this section we derive the relation r = A/(d(a - 1)), first given in 0. First, we define the phenotype robustness r 
as a weighted average over genotype robustness; second, we define it as the extent to which mutation off the neutral 
network does not deplete the growth due to fitness. Readers who are less interested in the biological motivation may 
skip this section. 

Genotype and phenotype robustness. Consider a neutral network P, and let its adjacency matrix be A. The genotype 
robustness r,- [T| of a genotype g, is the probability of a mutation being neutral: the number of neutral edges incident 
to gi (i.e. edges which do not lead to a different phenotype) divided by the total number of incident edges d(a - 1). 
The genotype robustness can therefore be written as 


_ ^jAij 
' d(a - 1) 


( 1 ) 


For a neutral network, let n(t) be its population vector at time f, with the /th component n,(t) corresponding to the 
population on genotype g,. The normalized population is distributed according to n(t)/ E; Suppose for now that 
in the limit f — » oo, the normalized population is distributed according to a unique distribution. Then we define the 
phenotype robustness r to be the long-time population-weighted average of the genotype robustnesses r,: 


E/«/(°°) " 


( 2 ) 


It is the fraction of the population flux that is neutral. We will now determine this limit. 

Mutational flux and fitness. Mutation induces a population flux across neighbouring genotypes. If the mutation rate 
per letter is p, the mutational flux is 1 - (1 - p) d(a 11 ~ pd(a - 1) for pd(a - 1) <s 1. It is the fraction of a population 
that mutates per generation. Some of this mutational flux will also cross phenotypic boundaries when neighbouring 
genotypes lie in two different phenotypes. That which does not cross phenotypic boundaries is neutral. The fitness / 
is the raw reproductive rate of the phenotype. After t generations, the total population of a neutral network will have 
changed by a factor of /', in the absence of mutations. 

Mutation matrix. The action of mutation on the population distribution over a single generation can be expressed 
by the mutation matrix M: 


M = (1 - pd(a - 1))/ + pA. (3) 

The first term is the probability that no mutation occurs and the second the probability of mutating. Being symmetric, 
A can be diagonalised by an orthonormal set of eigenvectors xg. 

M = (1 - pd(a - 1)) ^ XjxJ +p^j x i y i x J’ (4) 

i i 

where the x, satisfy the eigenvalue equation Ax, = v,x,. The population vector n(t) is obtained by transforming an 
initial vector no by M‘ and multiplying it by f: 

n(t ) = f'M’n(O) = 2 xjn(0)f (l - pd(a - 1) (l - JJ x ( . (5) 
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Let vi be the largest (principal) eigenvalue of A, denoted hereafter A. Since |v,j < d, all terms i > 1 decay exponentially 
with respect to the first for /r > 0. In the large time limit the sum is dominated by the first term, whose eigenvalue 
vi = A is largest: 

n, « xjn(0)f‘ |l -pd(a - 1)^1 - | x,. (6) 

We now show that defining the robustness as A/(d(a - 1)) agrees with the definition of phenotype robustness in (J2|. 
Indeed, plugging ([T]) into & 

2,;A, 7 H/(oo) 

r = - 

d(a - 1 )Z,«,(°°) 

YjjAn j(oo ) 

d(a - 1 )2i«i(°o) 

A 

d(a - 1) 

The quantity r therefore measures how well the shape of the neutral network can reduce the rate of deleterious 
mutation acting on the population as a whole. We see from 0 that at large time t at every generation, a fraction 
pd(a — 1)(1 — r) of the population mutates off the neutral network, and the growth rate (1 - pd( 1 - r))f is the fitness 
that can be usefully employed to increase the population and not spent replenishing population lost to deleterious 
mutations incurred at the boundary. The steady state distribution of the population depends only on the shape of the 
neutral network and on neither the mutation rate p nor the fitness /. 

2. Neutral networks with large eigenvalues 

Bricklayer’s graphs. Just how robust a phenotype can be—or how large an eigenvalue a neutral network can have— 
has remained an open question. For short sequences (d < 4, a — 2), we found from exhaustive enumeration that 
the most robust neutral networks are themselves hypercubes or interpolations between them, illustrated in Figure 1. 
Computational sampling for slightly longer sequences (5 < d < 9, a — 2) agrees with this. We generalize the sequence 
of graphs and interpolations between them in Figure 1 for a > 2 as follows: suppose all vertices { q} in H,i „ are labelled 
as integers from 0 to a d - 1, and two vertices share an edge if their base a representations differ in exactly one digit. 
Then G, u , is the subgraph induced by the vertices q < n. We call these graphs G, ul “bricklayer’s graphs” because they 
form the sequence by which a bricklayer would instinctively fill in the Hamming graph H r j a . For the remainder of this 
Note we set the alphabet size a - 2, so we are only concerned with hypercubes and their subgraphs. For simplicity 
we denote G , h2 by G n . We conjecture an extension of our main result for general a in the Conclusion. 

For neutral networks on sequences of short length d, the bricklayer’s graphs G„ are the most robust; they have 
maximal principal eigenvalues. For longer lengths d , however, a surprise is in store: the G„ are not the most robust 
neutral networks. In particular, we discovered the following counterexamples for d — 19 and above. Let S n be the 
star graph: a tree with one internal vertex and n leaf vertices. The principal eigenvalue of S„ is readily found to be 
sfn. Now let us compare the eigenvalue of a star of n vertices, S'i, to the eigenvalue A n of G„. As we prove below, 
A n < log 2 n. For the bricklayer’s graphs G„ to win, we need V/; - 1 < A n , implying 

- 1 < log 2 n. 

However, this not true for n > 20. It is an open question as to what shape does maximize the eigenvalue for larger 
graphs. In the concluding remarks of 0, the authors consider the possibility that Hamming balls (graphs consisting 
of all points that are most a given distance from a point) are asymptotic maximizers of hypercube subgraphs, but then 
provide some evidence that they are not. 

Our main result. In this Note we prove that the principal eigenvalue A„ of the bricklayer’s graph G„ satisfies 
A„ < log 2 n. Our general approach is to show by a geometric staircase argument that for d > 3, a slightly stronger 
inequality ( A„ < log 2 (n - 1)) holds for most n; it will then suffice to examine the cases where n - 2 d + 1, using 
polynomials that have A^d ± \ as roots. 
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Theorem. For all graphs G„, we have A n < log 2 n, with equality if and only if n is a power of 2. 


Equality is attained in the Theorem if n is a power of 2 since A„ must lie between the mean and maximum vertex 
degree 0 , and for n a power of 2, all vertices are of degree log, n. We wish to show that if n is not a power of 2, there 
is strict inequality. We make two observations. Observation 1: Since the principal eigenvalue of a proper subgraph 
of a connected graph is less than the principal eigenvalue of the graph itself, it follows that if n < m then A„ < A m . 
Observation 2: Since Gin = G„OK 2 , and the spectrum of a Cartesian product of graphs is the sum of their individual 
spectra 0 , it follows that if A„ < log 2 n then A 2n = A„ + 1 < log 2 2 n. Using these observations, we claim: 

Lemma 2.1. The Theorem is true for all n if for some k, 

A„ < log 2 0z - 1) for 2 k + 2 < n < 2 k+l - 1, (7) 

and also 

A 2 d^\ < \og 2 (2 d - 2), d > 5, and (8) 

A 2 d + 1 < log 2 (2 d + i), d > 3. (9) 

Proof The “staircase” argument is illustrated in Figure[2] We verify numerically that the Theorem is true for n < 16 
and 0 holds for k = 3. Now if 0 is true for some k, then by Observation 2, 

A 2n < log 2 (2 n - 2) for 2 k + 2 < n < 2 k+l - 1. (10) 

By Observation 1 we have the expansion 

hn-2 < 'hn -1 < Mn < log 2 (2n - 2) < log 2 (2« - 1) < log 2 2n (11) 

for 2 k + 2 < n < 2 k+ 1 - 1, so that A m < log 2 m for 2 k+i +2 < m < 2 k+1 - 2. Conditions 0 and 0 then ensure 

that A m < log 2 m for m = 2 k+ 1 + 1 and m = 2 k+2 - 1 as well, proving the Theorem for 2 k+ ~< n < 2 k+2 . Finally, 
note that CD together with 0 and 0 imply that 0 holds with k replaced by k + 1, so we may repeat our induction 
indefinitely. □ 

Therefore, the Theorem reduces to 0 and 0. which we will prove by looking at polynomials that have the 
eigenvalues of our desired graphs as roots. We will make use of the following standard theorem in linear algebra: 

(Cauchy’s Interlacing Theorem Q). Let A be an n X n symmetric nonnegative matrix with eigenvalues a \ < ■ ■ ■ < a n , 
and let B be an mxm principal submatrix of A with eigenvalues b\ < ■ ■ ■ < b m . Then for all j < m+ 1, Uj < bj < a„- m+ j. 


3. Polynomials with eigenvalues d 2 ^ ± i as roots 

Let Xn be the characteristic polynomial of the adjacency matrix of G„. Enumerating the hypercube spectrum, we find 


X 2 *(x) = Y\(x-(d- 2/))®, 


/=0 


so we may define the polynomial 

P 2 d(x) 


X 2 «(x) 


nt!( x -(d-2i)p- 


7\{ X -(d-2i)). 


( 12 ) 


(=0 


By applying Cauchy’s Interlacing Theorem to the adjacency matrices of G 2 ,/ +1 and Gv (thus “sandwiching” the 
spectrum of xi d +i by that of^ 2 d), we see from the multiplicity of the eigenvalues (of^f 2 <i) that X 2 d +i( x ) has the factor 
nfc, 1 (x-(d- 2i)P ~ 1 . A similar argument shows that^' 2 rf_ 1 (jc) has the same factor. So we may define the polynomials 


P 2 d±i(x) = 


X2“±l(x) 


nt t '(x-(d-2i)p- 
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Furthermore, we see that^ 2 rf - i has at most d simple roots, all of which are also roots of the ^-degree polynomial Py'-i ■ 
The reason we use the polynomials P is related to the fact that Gyi_\ is a Hamming ball of radius d— 1 in !l,/y. We 
define B lkr , the ^-dimensional ball of radius r, as the set of points in I l,jy that are Hamming distance at most r from 
the “origin” (the point labelled “0” according to the labelling scheme specified in the first paragraph of section 2). We 
determine recursive equations that give A, the principal eigenvalue of B r ] t . Consider the corresponding eigenvector w, 
and note that by symmetry the component of w corresponding to a given vertex depends only on the distance of the 
vertex from the origin. Therefore, let w k be the value of the component of w corresponding to a vertex of distance k 
from the origin. By matrix multiplication, we find that 


£ 

II 

o 

£ 

(13) 

Aw k = kw k -1 + (d — k)wk-ir\ for 1 < k < r 

(14) 

Aw r = rw r - 1 . 

(15) 


By setting wq — 1 and following the equations above for each fixed r, we find that the principal eigenvalue of By r is a 
root of the polynomial p r (A ), where po = A, p i = A 2 - d, and 


p r = Ap r _\ - rid - r + l)p r -2 for r > 2. 


Applying this to we can generate polynomials in A with coefficients in d (say fk(d, /!)) such that when k is 
substituted for d, the resulting polynomial in A has Ay \ as a root. Then f\{d. A) = A, fiid. A) = A 2 —d and, in general, 

fk(d. A) = Af k . t (d, A) - (k- 1 )(d-k + 2 )f k . 2 (d. A), (16) 

and f k (k, A) has Ay \ as a root. In fact, f k (k, A) has every simple eigenvalue of GA-i as a root (by the reasoning of the 
derivation). Since the degree of f k (k,A) as a polynomial in A is d, and we have from above that xv-\ h as at most d 
simple roots (all of which are also roots of the c/-degree polynomial fV-iX it must be the case that 


Py-\(A)=f k (k,X). 


(17) 


n 

1 2 3 4 5 6 7 8 9 10 12 14 16 



Figure 2: The “staircase” argument for the proof of Lemma Because each bricklayer’s graph is a subgraph of its successor, knowing the 
principal eigenvalue of a bricklayer’s graph immediately places restrictions on higher-dimensional bricklayer’s graphs (as demonstrated by the 
staircase figures of each color). 
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4. Bounding A 2 d_ 1 

Rewriting the right side of ([8]i by applying the Taylor expansion with Lagrange remainder gives 


l° g 2(2 —2) = d + log 2 \ l - —j\> d + 


log 2 \ 2 d 2 2d - 1 )' 


l°g 2 (2 £/ - 2) > d - l 


( 18 ) 


Therefore, for d > 5, 

Now we deal with the left side of ([8]). 

Lemma 4.1. A 2 d_ 1 <d — P 2 ,i_ \(d)/Pf_fd). 

Proof. Note that the function P 2 d_fx) is convex on x > A 2 <i. \ ■ To see this, let fix) = f]"_| (x - r,) be any monic 
polynomial with all real roots and observe that 


/"(*) = 2 ^ p| (x- n) 


h Jz=l \-i=hi±h,i±h 
ji<h 

which is always nonnegative if x is at least the largest root of /. Now since the tangent linear approximation of a 
convex function is an underestimate, we obtain P 2 d_fd) + If , fd) ■ (A 2 d_i - d) < 0, which implies the lemma. □ 

We evaluate the desired values of P 2 <‘-\ and its derivative using (17 i and the recursive relation in {Ill- 

Lemma 4.2. P 2 d_i(d) = d\ 

Proof. Recall the definition of fifd. A) from the previous section and the fact that Px fA) = ffk. A). We seek to prove 
that fk(k, k) — k\, and to do this we will prove a stronger claim that for i, k e N + , 

fk(f i) = (f)k> 

where (i)h = i(i - 1 )(i - 2) •••(/- k + 1) is the Pochhammer symbol. We use induction on k. For k - 1 we have 
fi(i, i) = i and f 2 (i, i ) = r - i = i(i - 1). Supposing the claim is true for f- 2 (i, i) and fk-i(i, i), we find from (161 that 

fk(i, i ) = ifk-\(f i ) -(k- 1)0' -k + Z)f k - 2 (i, i) 

= i(i)k-\ ~(k~ 1)0 ~k + 2)0)a- 2 
= 0 )*. 


n 


Lemma 4.3. 


d-i 


P’ 2d _M) = d'Y J - 1 


V 


Proof. With respect to the polynomials f k (d, T), let f k {d , d) denote d). Then Pf-fA) = f k {k , d), and f[{d, d) = 

1, f 2 (d, d) = 2d, and for k > 3, from (16 1 , 

f'(d. A) = f k .fd, d) + Afifd, A) -(k- l)(d-k + 2 )f'_ 2 id. A). 

We seek to prove that f'Ak. k) = k\ 2to tl’ and to do this we will prove a stronger claim that for integers k > 0, i > 0, 


k -1 


f k (k + i,k + i ) = k\ ^ 

i=o 


k-j + i- l\ 2 J 

' . )jTT- 
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We use induction on k. For the base cases k — 1 and k — 2 we find that /'(l + /, 1 +/) = 1 and / 2 (2+/, 2+i) = 2(2+/), 
as desired. Now supposing the claim is true for f' k _ 2 {k + i,k + i ) and ^ (k + i,k + /), it follows that 


f[(k + i,k + i) = fk-i(k + i, k + i) + (k + i)f[_ t (k + i,k + i) - (k - 1)(7 + 2 )f k _ 2 (k + i,k + i) 

= (* + »Vi + (* + «■)(*- D! Z - (*- D! Z 

.,=0 J=0 

*r-2 k-2 

- S±fi[ . ft n (t+01 . V [ (*-!)! (t-M! j . ; (*-!)! (t-j+i-l)! ] jft _ /> _ ivt t . »+-!)! _ V Jftftlft ifti fcMl 

— <i+l)l (,•+!)! ^ / j \ (fc— 7 + 2 )! (i+ 1 )! -TJ(k-j- 2 )! (i+1)! J j+1 ^>(i+ 1)! Z_j (k-j- 2)! 7 (i+l)l 

/=0 7=2 

(*—l>! 2* I , (*-!)! 2+ 1 ] (*~2+Q! 

(*-7+2)! 7+1 T v j J (7+1)! 

=‘sf + i: [</a (ft - + 0 - t] w + w» -w-w-ws 

2=2 
k -2 

_ 7.2+0! V * ! (2ft _ Zft) (*-j+Q! , 1 / 2 t - 1 _ 2*- 2 ) 

(7+1)! Zj (*-2-1)! V2+1 j ) (7+1)! " rK ‘V r *-l / 


= £§ + (*-!)£§ + (*" 1)^-2)^# + 


2=2 

k -2 


=* ! Z[ 


7=0 

r-i 


2=2 

(*-7+0! 21 _ (*—2-l)(*—7+7-1)! 21 ~ 

(i+l)!(fc-7-l)! 7+1 (i+l)!(r-7-l)! 7+1. 


= ft! y <ftft±! 

■ Zj '!(*-/ 


-7+7+1)! 22 


7=0 


1-7+1)! 2+1 ' 


□ 


Now substituting the results of Lemmas 4.2 and 4.3 into Lemma 4.1 

i ft d — 


yd-l 2-i * 
+J 7+0 j+i 


(19) 


A simple induction with base cases 1 < d < 3 shows that the sum in the denominator of the second term of ( p~9] > 
the right-hand side can be bounded by 


on 


A 2j „ 2 d 

/ — < 3 —. 

^ 7 ^ 


1=1 ■' 


Therefore, 


1 2d 

— < d - j. 


X'd-l 2i T. Id' 

+22=0 j+1 

Since we may combine the bounds (21 1 and (18 1 for d > 5, we have proved (|8j. 


( 20 ) 


( 21 ) 


5. Bounding A 21+1 

Rewriting the right side of <[9]> in much the same way as the previous section, we find that for d > 3, 

1 

log 2 (2 d +\)>d+^ r (22) 

As for the left side of (|9]>, we have a computational shortcut: 
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Lemma 5.1. Let G be a graph, and let e be a bridge ofG (i.e., e is an edge such that removing it would increase the 
number of connected components ofG). Let G* be the graph G with e removed, and G** be the graph G with e and its 
endpoints removed. Then \c = Xc ~XG"> where Xg is the characteristic polynomial of the adjacency matrix ofG. 

Proof This is an extension of Lemma 1 in J8), which states the result when G is a forest and e is any edge, and is also 
Theorem 1.3 in J9|. The proof involves expanding the matrix whose determinant is^ c; , using Laplacian expansion and 
linearity of the determinant. □ 

Corollary 5.2. X2 d +i(A) = Axv(A) -X2 d -i(A), an d so dividing by nf=i '(x — (d - 2 i))© -1 , 

p* + lU) = AP 2 4A) - Py-M). 

Lemma 5.3. A 2 j+ i < d - Py +l (d)/Pf +l (d). 

Proof From the preceding corollary. 


|4T 

Tot 


Pf(x) > 0 for 


Py + 1 W = 2 P' 2d (x) + xP”,(x) - P' 2d _ x (x). 

We wish to show that this is nonnegative on x > d. From the argument of the proof of Lemma 
x > d, and from the equation for the second derivative of a polynomial there, P'f(x) > P" d _ (x) for x > d since the 
roots of P 2 d_ | interlace those of P 2 d by Cauchy’s Interlacing Theorem. So P 2 d +1 (x) is convex on x > d, and so by 
linear approximation, P 2 d +1 (d) + P' 2d+l (d) ■ ( A(G 2 d +1 ) - d) < 0, which implies the lemma. □ 

Now as in the previous section, we evaluate the desired values of P 2 d +l (d) and its derivative. 

Lemma 5.4. P 2 d+i(d ) = -d\. 

Proof Using Corollary |5.2| and Lemma 4.2 P 2 d +X (d) = dp 2 d(d) - p 2 d_ x (d) = 0 — d\. 

Lemma 5.5. 


□ 


P' 2 d +X (d) = d\ 


d 1 2 j ^ 


d2 d - V —- 

U J+l 


Proof. Differentiating (12 1 and then plugging in d gives P' nd {d) = 2 d d\. Then using Corollary |5.2| and Lemma 

P 2d+X (d) = P 2 d(d) + dP' v (d) - P' 2d _ x (d) 


4.3 


= 0 + d2 a d\ 


= d\ 


d 2 d 


d A 2 j 

d\ > 

U > +1 

y 2> ^ 


□ 


Substituting the results of Lemmas 5.4 and 5.5 into Lemma 5.3 


A 2 d + 1 < d + 


1 


d 2 d ■ 


2-^ ' 

^j =0 7+1 


(23) 


From (|20|, we obtain that for d > 1, 


d + 


< d + 


dV-Yf-Jtjk (d-±) 2“ 

Since we may combine the bounds ( [24] ) and ( [22] ) for d > 5, we have proved & 
This concludes the proof of the Theorem. 


( 24 ) 
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6. Conclusion 


As an additional remark, if we bound log 2 (2' / 
result: for any N, there exists D so large that A 2 j 


N ) for general N in the manner of (18 1 , we can deduce an asymptotic 
< log 2 (2 J 


i iug 2 r^ - N) for all d > D. Similarly generalizing (22 1 leads to 
the result that for any e > 0, there exists D so large that Ay+i < log, (2^ + e) for all d > D. 

Throughout most of this paper, we have set a - 2. We conjecture an extension of the Theorem for general a: 

Conjecture. For all graphs G„ M , we have A n , a < (a — 1) log n n, with equality if and only if n is a power of a. 


There are several interesting and potentially important questions that we have not considered here, which merit 
further investigation. We prove that the form of the hypercube subgraph with maximal eigenvalue is a bricklayer’s 
graph for small d but the general form of the maximizers is unknown. Indeed it is an open avenue of study to find 
even non-trivial bounds on the eigenvalue of a hypercube subgraph in terms of its number of vertices. While for 
small dimension d the bricklayer’s graphs are optimal, for ,v > 20 and d > s — 1, Hamming balls of radius 1 are 
superior. For large s and d <s s, Hamming balls of larger radius may eventually dominate, but this is unproven. How 
these transitions extend to larger values of alphabet size a is also an open question, though it seems that the critical 
dimension separating bricklayer’s graphs and balls grows with a. We hope that further research by others will shed 
light on these questions. 
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